Stationary light in cold atomic gases 
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We discuss stationary light created by a pair of counter-propagating control fields in A-type atomic 
gases with electromagnetically induced transparency for the case of negligible Doppler broadening. 
In this case the secular approximation used in the discussion of stationary light in hot vapors is no 
longer valid. We discuss the quality of the effective light-trapping system and show that in contrast 
to previous claims it is finite even for vanishing ground-state dephasing. The dynamics of the photon 
loss is in general non exponential and can be faster or slower than in hot gases. 

PACS numbers: 42.50.Gy, 32.80.Qk 



I. INTRODUCTION 

Strong coupling between light and matter is of large in- 
terest in many fields of physics. It is of particular impor- 
tance in quantum information and quantum-optical real- 
izations of strongly interacting many-body systems. The 
interaction strength between single photons and quan- 
tum dipole oscillators is determined by the value of the 
electromagnetic field at the position of the oscillator and 
thus by the spatial confinement of the photons. This 
has lead to the development of cavity quantum electro- 
dynamics where strong confinement and thus strong cou- 
pling is achieved by means of low-loss micro-resonators 
PI . An alternative suggested by Andre and Lukin 0] and 
first implemented in a proof-of-principle experiment by 
Bajcsy et al. ]3f, is to create spatially confined quasi- 
stationary pulses of light with very low losses by means 
of electromagnetically induced transparency (EIT) 0, Q 
in an ensemble of A-type three-level atoms driven by two 
counter-propagating control fields. The physical proper- 
ties of stationary light were discussed in a number of the- 
oretical studies. It was shown that under adiabatic con- 
ditions quasi-stationary light obeys a Schrodinger equa- 
tion with complex mass and that inhomogencous control 
fields can be used to spatially confine and compress its 
wave-packet @. The fundamental quasi-particles of sta- 
tionary light have been identified [7] , a transition from a 
Schrodinger-like to a Dirac-like dynamics has been found 
Q and many-body phenomena with stationary-light po- 
laritons have been discussed theoretically such as the 
Tonks gas transition Q and Bose-Einstein condensation 

An essential assumption of the original model for sta- 
tionary light is the secular approximation in which spatial 
modulations of the ground-state coherence of the A-type 
atoms with wavenumbers on the order of the optical fields 
and its harmonics are neglected. The latter is a very good 
approximation in warm gases, where atomic motion leads 
to a fast dephasing of fast spatial oscillations. It fails 
however for cold gases or other systems where the mo- 
tion is suppressed such as solids or atoms in optical 
lattices [13, [H| . Although the problem of a secular ap- 



proximation can be entirelyavoided by using a double-A 
rather than a A transition 0, HH > it is interesting to con- 
sider the dynamics in a cold gas of A-type atoms without 
the secular approximation. The earliest analysis of this 
case was done by Moiseev and Ham 15]. However their 



analysis was limited to the case of unequal coupling field 
intensities, thus the probe field was not really station- 
ary. In a more recent study M0lmcr and Hansen found 
that without secular approximation the wave-packet of 
light is truly stationary for equal strength of the control 
field, i.e. does not undergo a diffusive spreading, and 
the only source of photon loss is the finite lifetime of the 
ground-state coherence [161 ] . In this analysis radiative 
losses where neglected however. The result obtained in 
fl(| predicts the possibility of light trapping in EIT me- 
dia with unprecendented Q factors. In the present paper 
we analyze stationary light in A-type media without the 
secular approximation by taking into account the relax- 
ation of the excited state. We prove that the unavoidable 
excited-state decay limits the lifetime of the probe field 
in the medium. It leads either to a broadening of the 
quasi-stationary wavepacket in time or a splitting into 
two pulses fPi\ depending on the system parameters. The 
general dynamical behaviour is non-trivial, leading e.g. 
to a non-exponential decay of photons from the initial 
volume. We identify parameter regimes in which the ef- 
fective loss in cold gases is slower or faster than the one 
in a warm gas where the secular approximation holds. 



II. FIELD EQUATIONS OF STATIONARY 
LIGHT BEYOND THE SECULAR 
APPROXIMATION 

We here consider a medium consisting of an ensemble 
of non-moving three-level atoms with a A configuration 
shown in Fig[T] We assume that initially some coherence 
is stored in the lower levels of the atomic medium, so 
that when a standing wave resonant coupling field f2 c 
is applied, a quasi-stationary probe field E is created. 
For simplicity the states \g) and \s) are assumed to be 
degenerate, thus the wave vectors of the probe and the 
coupling fields have equal magnitude k. 
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FIG. 1: (color online) Schematic diagram of the three-level 
atomic system and the medium pumped by a standing wave 
coupling field. 



The interaction Hamiltonian in rotating wave approx- 
imation reads 

H = -hY,(gE<Ti g + n c ai s ) + H.c, (1) 
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where E is the dimensionless slowly- varying complex am- 
plitude of the probe field, g is the atom field coupling 
constant, fi c is the Rabi frequency of the coupling field, 
and ct^jj — \fJ>)jj {v\ are the atomic transition operators of 
the jth atom between states and \i>). The atom dy- 
namics is governed by Langevin equations corresponding 
to U) and including losses from the excited state. They 
can be written as 

dtfrge = -r<7 ge + igE(a gg - <7 ee ) + iO, c a gs , (2) 
d t a gs = in,*a ge + igEa es , (3) 

where T is the relaxation rate of the upper level and it 
is assumed that the decoherence of the lower-level tran- 
sition is negligible on the time scale of interest. In the 
limit of low probe-field intensities (\gE\ <C |^ c |) an d for 
an initial preparation of the ensemble in the ground state 
\g), we can set in eq.© a gg w 1; a ss = a ee = a se = 0, 
which corresponds to the well-known pertubative linear- 
response limit. Since losses are included in the above 
equations there should be in general Langevin noise oper- 
ators [l8j . The noise operators are however inconsequen- 
tial since in the linear-response limit, considered here, 
there is no excitation of the excited states. Thus they 
are neglected. We assume furthermore, that the charac- 
teristic duration of interaction T is long compared with 
respect to the upper level relaxation (FT 3> 1). This al- 
lows for an adiabatic elimination of the optical coherence 



a ge and equations ([2]) are reduced to 

TlTge = igE + i$l c Vgs (4) 

d t a gs = ifll&ge- 

Differentiating the first equation with respect to time and 
assuming a constant control field f2 c yields 

OtCTge = Y"t E — a gei 

which has the formal solution 

= jLJ _exp|-^(t-T)|dr 

o 

+«^fo- ffS (^O)exp | -l^-tj . (5) 

Since the coupling field is a standing wave formed by two 
counterpropagating fields of equal intensity and polariza- 
tion, it can be expressed as £l c (z) = 2f2cos(fcz), where 
Q represents the amplitude of the coupling field. The 
probe field consists also of two counterpropagating com- 
ponents E — E + e lkz + E_e~ lkz . Due to the presence of 
|fi c (z)| 2 in the exponents in eqj5l the optical coherence 
and thus by virtue of eq|4] also the ground-state coher- 
ence will develop all harmonics of e ±lkz . Thus we make 
the ansatz 

oc 

a gs (z,t) = J2 *i?(z,t)e iknZ , (6) 

n— — oo 
oo 

a ge (z,t) = J2 

n— — oo 

The secular approximation corresponds to disregarding 
all terms in a gs with n ^ 0. This is justified in a hot gas 
where atomic motion washes out the fast spatial oscilla- 
tions associated with terms and n ^ 0] i @ • 

The propagation of the probe pulse components are 
governed by the Maxwell equations for the slowly varying 
field amplitudes 

where cx^ 1 ^ are the components of atomic coherence be- 
tween levels \g) and |e) that oscillate in space according 
to e ±tkz , and N is the number density of atoms. 

If the stationary light pulse is generated from a stored 
spin coherence without rapidly oscillating components, 
i.e. for agt n \z, t = 0) = 0, for u ^ the corresponding 
initial conditions are 

Td£(z,0) = ina gs (z,0), (8) 
a±"(z,0) = 0, for n ^ 1. 
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Using the identity e Acosz = J2 I\n\ (A) cos (zn) = 

n— — oo 

oo 

I\ n \ (A) e tzn , with I n being the nth order modified 

n— — oo 

Bessel function, we can rewrite equation ([5|). 



(9) 



dE. 



hiP{t-*))-Qjr-hm-H)) 



+ia gs (z, 0)Cle 



-fit 



hlfit) - hi fit 



dE T 
~dt r 

/r, 



where we have introduced ft = 2£1 2 /Y. If we consider 
times t which are sufficiently large, the initial value term 
in eq.© can be disregarded. 

Substituting © and ([S]) into equations ([7]) and in- 
troducing the sum and difference normal modes E s — 
E++E-, E d = E+- E- yields 



dt dz 



dE d 
dt 



dE. 
dz 



g 2 N 



g^N 
Y 



dt 



(10) 

(ii) 



where f±(x) = e x Iq\x) ± h(x) . 

For the following discussion it is convenient to intro- 
duce normalized variables and parameter 



9 2 JV 



fc abs 



2fl 2 
9 2 N 



2 cot 2 6, (12) 



where Z a b s = cY/g 2 N is the resonant absorption length 
of the medium in the absence of EIT. This leads to the 
normalized equations 



d T E s + d^E d = - / dr' f_(a(T-T r ))d T ,E s , (13) 
Jo 

d T E d + 8^E S = - [ dr' f+(a{ T - r')) d T ,E d . (14) 
Jo 

III. PROBE-FIELD DYNAMICS 



In the following we will qualitatively discuss the probe- 
field dynamics resulting from eas. (|13I14|) . illustrate the 
results with numerical examples and compare the field 
evolution with the case of a hot atomic gas. Eqs. (|13ll4p 
turn into the corresponding equations for a warm atomic 
gas where the secular approximation is valid, if one sets 
f-(a(r - t')) ~ (2/a)<S(r - r'), and f+(a(r - r')) - 1 



dE s 
dr 
dE d 



dE d 



d£ 
8E S 
dr ' «9£ 



tan 2 9 



-Eh 



dE s 
dr ' 



(15) 
(16) 



where tan 2 9 = g 2 N/tt 2 . In this case adiabatic eliminat- 
ing the fast decaying difference mode, i.e. E d ~ — d^E s , 
results in a diffusion equations for the sum mode 



or 



dE s 
dr 
dE s 
dt 



d 2 E s 



0. 



o 2 e s 

dz 2 



= 0, 



(17) 
(18) 



where u gr = ccos 2 is the group velocity of EIT. As- 
sociated with the diffusion is a (non-exponential) loss 
of excitation with a characteristic time scale of T* oss = 
L 2 /Zabsfgr, with L being the characteristic initial confine- 
ment length of the stationary pulse. 

In order to discuss the stationary-light dynamics be- 
yond the secular approximation we start with numerical 
solutions of eqs. (|13U4[) for two characteristic cases. In 
Figf2]the decay of the total field intensity I = j dz(E 2 _ + 
i? 2 ) in the interval {— 3Lo,3Lo} is shown after retrieval 
of an initial gaussian spin excitation of spatial shape 
exp{— z 2 /Lg}, and Lq — 5/ a b s (solid line) for two im- 
portant cases. In the first case (top curve) tan 2 9 = 100, 
i.e. a = 0.02, in the second (bottom curve) tan 2 9 = 1, 
i.e. a — 2. Also shown is a comparision with the results 
obtained with the secular approximation (dotted line). 
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FIG. 2: (color online) Decay of total field intensity within the 
spatial region from -3 Lo to 3 Lo for an initial spin excita- 
tion of spatial shape exp{— z 2 /Lq}, Lo = 5i a b s in cold atomic 
gas (blue, solid line) and in inhomogeneously broadened (red, 



dotted line) media for tan 2 9 = 
(bottom). 



100 (top) and tan 2 6 = 1 



The time evolution of the field distributions of E s and 
E d for the two cases are shown in Fig[3] (tan 2 9 = 100) 
and FigH (tan 2 9 = 1). 

From the numerical examples several conclusions can 
be drawn: First of all one recognizes that contrary to 
the claims in Ref.[l6[ the field intensity decays even if 
the dephasing of the ground-state coherence is neglected. 
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FIG. 3: (color online) Spatio-temporal dynamics of sum E s 
(a) and difference Ed x f0 (b) modes of the probe field corre- 
sponding to the top curve in Fig[2j i.e. tan 2 9 = 100. 



The decay is caused by the relaxation of the upper state 
which was not taken into account in (lfl | by restricting 
the discussion to the loweest order in the adiabatic ex- 
pansion. Thus stationary light in cold gases or solids 
does not provide a perfect cavity. Secondly the decay of 
the intensity can either be slower or faster as compared 
to the case with secular approximation. In the first case, 
i.e. Fig|3] the evolution of the field distribution is very 
similar to the diffusive spreading but much slower than 
in the secular-approximation limit shown in FigO On 
the other hand in the second case (see FigSJ) , two pulse 
components emerge which propagate with the group ve- 
locity u gr = c cos 2 9 with some additional loss [17| • 

We now want to give a quantitative explanation of the 
different dynamics in the two cases, which is due to the 
different action of the integral kernels f± in eqs. ([T5)) and 
(Q3J). For this it is instructive to perform a Laplace- 
transform of cqs. (fT3|) and (I14[) : 



sE s - E S (Q) + d^E d = -/_(s) ■ [sE s - E s {0)) ,(19) 
sE d - E d (0) + d s E s = -f+(s) ■ [sE d - £ d (0)) (20) 

where the Laplace transform of f± reads 

1 
a 

1 

a 



/+(*) 
/-(«) 



2- + 1 - 1 
1 



This yields 



sE s -E s (0) 
sE d - E d (0) 



v/2f+T 

d^Ed 



1 



1 + /-(')' 
'! + /+(«)' 



(21) 
(22) 
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FIG. 4: (color online) Spatio-temporal dynamics of sum E 3 
(a) and difference Ed X 10 (b) modes of the probe field, as 
well as field component E+ (c), corresponding to the bottom 
curve in FigfJJ i.e. tan 2 — 1. The white lines indicate the 
velocity 2c cos 2 9/(2 cos 2 9 + sin 2 9). 



1. limit of small a 

Assuming that a is small, a series expansion of 1/(1 
f±(s)) yields 

1 s 



l + /±(«) 



1 



For the slow time evolution, i.e. for (physical) times large 
compared to (g 2 N/T)^ 1 only values of s < 1 are rele- 
vant, and thus the right hand side can be replaced by s. 
Substituting this into eqs. fTfj)) and (|2"01) one arrives at 



d 2 6 E s 
djE d 



E s = -(E 3 (^,0) + d s E d ^,0)) 
E d = -(E d (^0)+d i E s ^,0)) 



which describes truly stationary wave packets that decay 
exponentially with increasing distance. That there is no 
dynamics is of course due to the fact that only the leading 
order term in the expansion of 1/(1 + f±) was taken into 
account. 



2. limit of large a 

A qualitative explanation of the opposite case can be 
found by considering the limit of large a. To properly 
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FIG. 5: (color online) Dynamics of the sum E s (a) and dif- 
ference modes Ed X 10 (b) of the probe field in inhomogc- 
neosly broadened medium, i.e. under conditions that justify 
the secular approximation. All parameteres are identical to 
the Fig(2] 



In the limit of large a one arrives at wave equations 
for the forward and backward components 



d 2 T E±± — — zd 2 E± = 0, 

(a + 1 — 1 jx) (a — 1 + x) 5 



which reads in physical time and space: 



d 2 E ± ± 



(a + 1 - l/x)(a - 1 + x) 



d 2 E± = 0. 



Thus the envelope of E± evolves freely This explains 
the splitting of the stationary light wavepacket into two 
components each of which propagating with the mod- 
ified group velocity 2c cos 2 9/(2 cos 2 9 + / sin 2 9), with 
/ = (x — l/x)/4. Noting that the most relevant Laplace 
frequency for the example in FigH] leads to a value of / 
on the order of unity we find reasonable agreement with 
the numerical results. 



analyze this case one has to take into account that also 
the most relevant Laplace frequency s increases when a 
becomes large. In fact the numerical data suggest that 
the ratio a/s, with s being the most relevant Laplace 
frequency approaches a constant. Thus in this case one 
has 

1 a 



l + /+(s) 0-1 + a; 
1 a 



1 + /-00 a+l-l/z 



where x = \ 2 — hi is well approximated by a constant. 
This leads to the approximate equations 

, ,° u d 6 E d = 0, 
a + 1 — 1/x 

a 



d T E s 
d T E d 



a — 1 + x 



-dtE s 



0. 



IV. SUMMARY 



We considered the dynamics of stationary light in a 
standing medium without secular approximation and de- 
rived equations describing the evolution of the sum and 
difference modes of the pulse. A numerical as well as ap- 
proximate analytical solution showed that for small cou- 
pling field intensities the probe field spreading is slower 
than in the secular approximation but in contrast to the 
results of non-zero. In the opposite limit of strong 
coupling the probe pulse splits into two countcrpropagat- 
ing components. 

G.N. acknowledges support by the Alexander von 
Humboldt Foundation. 
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